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Abstract 

We calculate in chiral perturbation theory the double-pion photoproduction processes 

7r~7 — )• ir^TT^TT^ and -7r^7 — )■ ir^ir^ir^ at low energies. At leading order these reactions 

are governed by the chiral pion-pion interaction. The next-to-leading order corrections 

arise from pion-loop diagrams and chiral-invariant counterterms involving the low-energy 

O . constants ii, £2, ^3 and £4. The pertinent production amplitudes Ai and A2 depending 

on five kinematical variables are given in closed analytical form. We find that the total 

cross section for neutral pion-pair production 7r~7 — ?• 7r~TT^TT^ gets enhanced in the region 

^/s < 7m,7r by a factor 1.5 — 1.8 by the next-to- leading order corrections. In contrast to 

Q^ . this behavior the total cross section for charged pion-pair production 7r~7 — )• -k'^tt'-k' 

remains almost unchanged in the region ^/s < Qm.,^ in comparison to its tree-level result. 

. Although the dynamics of the pion-pair production reactions is much richer, this observed 

Oj. pattern can be understood from the different influence of the chiral corrections on the 

JL \ pion-pion final state interaction (tt+tt" -^ tt^tt'^ versus tt'tt" — )■ vt'tt"). We present 

Q ' also results for the complete set of two-pion invariant mass spectra. The predictions of 

i-^' chiral perturbation theory for the 7r~7 — )• Stt processes can be tested by the COMPASS 

experiment which uses Primakoff scattering of high-energy pions in the Coulomb field of 

^ a heavy nucleus to extract cross sections for 7r~7 reactions with various final states. 
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The pions (tt^, 7r°, tt^) are the Goldstone bosons of spontaneous chiral symmetry breaking in 
QCD: SU{2)l X SU{2)fj — )• SU{2)v- Their low-energy dynamics can therefore be calculated 
systematically (and accurately) with chiral perturbation theory in form of a loop-expansion 
based on an effective chiral Lagrangian. The accurate two-loop prediction [1] for the isospin- 
zero S-wave Trvr-scattering length Og = (0.220 ± 0.005)m~^ has been confirmed in the E865 [2] 
and NA48/2 [3] experiments by analyzing the vr+vr" invariant mass distribution of the rare kaon 
decay mode K^ — ^ 7i~^7i~e'^i'e. One particular implication of that good agreement between 
theory and experiment is that the quark condensate (0|gg|0) constitutes the dominant order 
parameter [4] of spontaneous chiral symmetry breaking (considering the two-flavor sector of 
QCD). Likewise, the DIRAC experiment [5] has been proposed to determine the difference of 
the isospin-zero and isospin-two S-wave Trvr-scattering lengths Oq — a^ by measuring the life time 
(r ~ 3fs) of pionium (i.e. tt+tt^ bound electromagnetically and decaying into tt^tc^). In the 
meantime the NA48/2 experiment [6] has accumulated very high statistics for the charged kaon 
decay modes K"^ — )■ n'^n^ii'^, which allowed to extract the value Og — Og = (0.268 ± 0.010)m~^ 
for the TTTT-scattering length difference from the cusp effect in the 7r°7r° mass spectrum at the 
TT+TT^ threshold. This experimental result is again in very good agreement with the two-loop 
prediction Og — Og = (0.265 ± 0.004)m~^ of chiral perturbation theory [1]. For a discussion of 
isospin breaking corrections which have to be included in a meaningful comparison between 



theory and experiment, see ref. [7]. Clearly, these remarkable confirmations give confidence 
that chiral perturbation theory is the correct framework to calculate reliably and accurately 
the strong interaction dynamics of the pions at low energies. 

Electromagnetic processes offer further possibilities to probe the internal structure of the 
pion. For example, pion Compton scattering 7r^7 — )■ 7r^7 at low energies allows one to extract 
the electric and magnetic polarizabilities (a^ and 13^^) of the charged pion. Chiral perturbation 
theory at two-loop order gives for the dominant pion polarizability difference the firm prediction 
a^ — P^ = (5.7 ± 1.0) ■ 10~^fm^ [8]. It is however in confiict with the existing experimental 
results from Serpukhov «„-/?„ = (15.6 ± 7.8) • 10"^ fm^ [9] and MAMI a^ - /3^ = (11.6 ± 
3.4) • 10~^fm^ [10] which amount to values more than twice as large. Certainly, these existing 
experimental determinations of aj^ — (3^^ raise doubts about their correctness since they violate 
the chiral low-energy theorem notably by a factor 2. The chiral low-energy theorem [11] 
relates a^r — /Sir = a(^6 — h) / {'^4:7i'^ f^rrij,) + 0{m.„) to the axial-vector-to-vector form factor 
ratio hA/hv = 0.443 ± 0.015 = (4 - 4)/6 + 0{ml) measured in the PIBETA experiment [12] 
via the radiative pion decay 7r+ — ?> e^u^'y. The two-loop calculations of refs.[8, 13, 14] assure 
that the 0{mj^) corrections to it are in fact small. 

In that contradictory situation, it is promising that the ongoing COMPASS experiment 
[15] at CERN aims at remeasuring the pion polarizabilities, Q;,r and Ptt, with high statistics 
using the Primakoff effect. The scattering of high-energy negative pions in the Coulomb field 
of a heavy nucleus (of charge Z) gives access to cross sections for n~'j reactions through the 
equivalent photon method: 

da _ Z^a Q^-Qi,n^ (. ^ _ s-ml 



Here, Q denotes the momentum transferred by the virtual photon to the heavy nucleus of charge 
Z, and one aims at isolating the Coulomb peak Q — )■ from the strong interaction background. 
The last factor 0-^-^(5) is the total cross section for a 7r^7 reaction induced by real photons with 
^/s the corresponding 7r~7 center-of-mass energy. Note that eq.(l) applies in the same form to 
differential cross sections on both sides. The COMPASS experiment is set up to detect simul- 
taneously various (multi-particle) hadronic final states which are produced in the Primakoff 
scattering process of high-energy pions. In addition to pion Compton scattering 7r^7 — > 7r^7 
(which is of primary interest for determining the pion polarizabilities a^r and (3^^) the reaction 
7r~7 — 7> 71^ Tc^ serves as a test of the QCD chiral anomaly (i.e. the anomalous V AAA rectangle 
quark diagram) by measuring the 737r coupling constant F^s^r = e/{4:7i^f^) = 9.72 GeV~^. For 
the two-body process 7r~7 — )■ 7r~'ir^ the one-loop [16, 17] and two-loop corrections [18] of chiral 
perturbation theory as well as QED radiative corrections [19] have been worked out. Thus an 
accurate theoretical framework is available to analyze the upcoming data. The 7r~7 reaction 
with three charged pions in the final state is used by the COMPASS collaboration in the energy 
range 1 GeV< y/s < 2.5 GeV to study the spectroscopy of non-strange meson resonances [20] 
(ai(1260), 02(1320), 712(1670), 7r(1800), a4(2040) etc.) and to search for so-called exotic me- 
son resonances [21] (e.g. 7ri(1600)) with quantum numbers different from simple (constituent) 
quark-antiquark bound states. The statistics of the COMPASS experiment is actually so high 
that the event rates with three pions in the final state can even be continued downward to the 
threshold. The cross sections (and other more exclusive observables) of the 7r~7 — )■ Stt reac- 
tions in the low-energy region ^/s < 1 GeV offer new possibilities to test the strong interaction 
dynamics of the pions as predicted by chiral perturbation theory. The total cross sections for 
the processes 7r~7 — > Tr^Tr^Tr" and 7r^7 — > tt+tt^tt^ at tree-level have been calculated recently 
in ref. [17] and it was found that the cross section for the charged channel (tt+tt^tt^) comes 



out about a factor of 2.5 larger than the one for the neutral channel (7r~7r°7r°). In both cases 
it is the chiral pion-pion interaction (together with the electromagnetic pion-photon coupling) 
which governs these reactions at leading order. A preliminary analysis [22] of the COMPASS 
data for 7r^7 — )■ n~^7i~7i~ in the low-energy range 0.5 GeV < a/s < 0.8 GeV leads to a total 
cross section which seems to be in agreement with the tree-level result [17] of chiral pertur- 
bation theory. An arguable aspect of that analysis is the absolute normalization which has to 
known well in order to convert count rates of events into a cross section. However, before any 
conclusions about an agreement between theory and experiment can be drawn the predictions 
of chiral perturbation theory need to be sharpened further by including higher orders in the 
small momentum expansion. 

This is precisely the purpose of the present paper. We evaluate for the double-pion photo- 
production processes 7r~7 — > tt~7[^tt^ and 7r~7 — )■ 7r'^7r~7r~ the next-to-leading order corrections 
as they arise from pion-loop diagrams and chiral-invariant counterterms (proportional to the 
low-energy constants ij). Our paper is organized as follows. In section 2, we treat first the 
somewhat simpler case of neutral pion-pair production. We give for individual diagrams the 
analytical expressions for the pertinent production amplitudes, Ai and A2, which depend on 
five independent Lorentz-invariant kinematical variables. In this detailed exposition one can 
check the exact cancellation of ultraviolet divergences between pion-loops and counterterms. 
Combined with the tree-level amplitudes, Ai and A2 are employed to calculate the total cross 
section and the two-pion mass spectra for 71^7 — )■ tt^tt^tt^. We find that the total cross sections 
(as well as the two-pion mass spectra) get enhanced by a factor 1.5 — 1.8 after including the 
next-to-leading order chiral corrections. Although the dynamics of the whole process is much 
richer this enhancement can be understood (in an approximate way) from the tt^tt" — )■ tt^tt^ 
final state interaction. In section 3, the same analytical calculation is performed for charged 
pion-pair production 7r^7 — )> tt+tt^tt^. In that case many more diagrams do contribute since 
the photon can now couple to all three outgoing (charged) pions. For that reason we restrict 
ourselves to specifying the finite parts of the pion-loop diagrams and the total counterterm 
contribution (rewritten in terms of the low-energy constants ij which subsume the chiral log- 
arithm ln(m^/A) generated by pion-loops). Interestingly, we find that the total cross section 
for 7r^7 — )> TT+TT^TT^ in the low-energy region SttItt < ^/s < 6m^ remains almost unchanged by 
the inclusion of the next-to-leading order corrections. Although the dynamics of the process is 
much richer the (known) weak infiuence of chiral corrections on the isospin-two tt^tt^ — ?> tt^tt^ 
scattering length Og [1] can provide an explanation for this feature. In this case the tt+tt^ mass 
spectrum reveals better certain dynamical details which get averaged out in the total cross sec- 
tion. We estimate also the uncertainties of the observables for the 7r~7 — )■ Stt reactions which 
are induced by the present errorbars of the low-energy constants ij, and find about ±5%. 

In summary, we have calculated the processes 7r"7 — ?■ ti^ii^tt^ and 71^7 — ?> tt+tt^tt^ at next- 
to-leading order in chiral perturbation theory. The predictions for the total cross sections and 
other more exclusive observables can be tested soon by the COMPASS experiment at CERN. 

2 Neutral pion-pair production 

In this section we treat the neutral pion-pair production process: tt~{pi) +'~f{k,e) — ?■ 7r~(p2) + 
'^^{Qi) + '^^{Q2)- We choose for the (transversal) real photon the Coulomb-gauge in the center- 
of-mass frame, which entails the conditions e ■ pi = e ■ k = 0. These conditions imply that 
all diagrams for which the photon couples to the in-coming pion vr~(pi) vanish identically. 
Furthermore, it is advantageous to parametrize the special- unitary matrix-field U in the chiral 






Figure 1: Tree-level diagrams for tt 7 — t- vr tt'^tt" and tt+tt tt 
pions. Only the left diagram contributes to 27r'^-production. 



Arrows indicate out-going 



Lagrangian £7^^ through an interpolating pion-field vf in the form U = J I — vr^//^ + ir ■ vf//^. 
It has the consequence that no jAii and 7671 contact vertices exist at leading order. Under 
these assumptions one is left with one single non- vanishing tree diagram shown in Fig. 1. Let 
us recall the expression for the total cross section for tt" 



7—7-71 71^71^ at tree-level: 
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Here, 



{Pi 



kY denotes the squared center-of-mass energy, a = e'^/An = 1/137.036 is the 
fine structure constant, and /tt = 92.4 MeV stands for the pion decay constant. In comparison 
to eq.(16) in ref. [17] we have included (some) isospin breaking effects by distinguishing the mass 
of the charged pion m^ = 139.570 MeV from the mass of the neutral pion m^^o = 134.977 MeV. 
Since it turns out that these effects are very small (see Fig. 7) we will perform the whole 
calculation in the limit of isospin symmetry. In order to present the next-to-leading order 
corrections from chiral loops and counterterms one has to start from the general form of the 
T- matrix, which reads (in Coulomb-gauge): 

2er 



T 



J TT 



e-qiAi + €■ g2^2 



(3) 



In this decomposition Ai and A2 are two dimensionless production amplitudes which depend 
on s = [pi + k)"^ and on four other independent (Lorentz-invariant) Mandelstam variables: 



Si = {p2 + qiY, S2 = {p2 + q2Y, ti = {qi-kf, t2 = {q2~kf 



(4) 



This set is very convenient for describing the permutation of the two identical neutral pions 
in the final state via (si <-)■ S2, ti ^ ^2)- For most diagrammatic contributions Ai = A2, but 
there are a few exceptions which require their separate listing and specification. 

With the form of the T- matrix defined in eq.(3) and the kinematical variables introduced 
in eq.(4) the tree-level amplitudes read: 
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(5) 



„ -tl-t2 

Note that T^*'''^'^) is written with physical parameters (/^, m^l 

parameters as given by the chiral Lagrangian. It requires an extra renormalization contribution 

(see eq.(26)) to properly account for this difference. 



and not with leading order 



2.1 Diagraminatic calculation 

In this subsection we present analytical expressions for the amplitudes Ai and A2 as they arise 
from the next-to-leading order chiral loops and counterterms. We go step by step through 
the whole set of contributing diagrams. The one- loop diagrams (I), (II), (III) shown in Fig. 2 
include an additional TTTr-interaction but leave the pion-photon coupling (of the tree diagram) 
unchanged. One finds from diagram (I): 
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with the abbreviation 
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for standard ultraviolet divergence in dimensional regularization. Note that ^ is always ac- 
companied by the chiral logarithm ln(m^/A). The complex- valued loop function J{s) has the 
form: 
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Similarly, one gets from diagram (II): 
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and the contribution from diagram (III) follows by interchanging the two neutral pions: 
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The one-loop diagrams (IV), (V), (VI) shown in Fig. 3 include also an additional TTTr- 
interaction but the photon couples now to a charged pion inside the loop. One finds from 
diagram (IV): 
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Figure 2: One-pion loop diagrams for vr 7 — t- vr tt'^tt^. Their combinatoric factor is 1/2. 






Figure 3: One-pion loop diagrams for vr 7 — )■ vr vr vr . Their combinatoric factor is 1. 



where the complex-valued loop function G{s) has the form: 
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Considerably more lengthy are the expressions which one obtains from evaluating diagram (V): 
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Figure 4: One-pion loop diagrams for vr 7 — t- vr tt'^tt^. Their combinatoric factor is 1/2. 
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In the process of evaluation one encounters loop integrals over three pion-propagators with 
tensors up to third rank (/^, lfj,y, l^lylp) in the numerator. After the pertinent tensor reduction 
all the occurring scalar loop functions can be expressed as linear combinations of J(. . .) and 
G{. . .) with rational coefficient functions. Note that the quality Ai = A2 does no more hold 
for diagram (V) since the four-momenta qi and q2 of the two neutral pions appear in the loop 
in a non-symmetrical way. The permutational symmetry between both 7r° is restored by the 
contribution from diagram (VI): 
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Diagram (VII) in Fig. 4 generates a (constant) vertex correction to the pion-photon coupling 
and therefore leads to the amplitudes: 
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which are proportional to the tree amplitudes written in eq.(5). The chiral six-pion vertex 
appearing in diagrams (VIII) and (IX) represents some challenge with respect to tackling the 
combinatorics involved. Altogether, one finds from diagram (VIII) the momentum- dependent 
amplitudes: 
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Figure 5: Diagrams for 7r^7 — )■ Tr^vr'^Tr'^ involving the pion wave function renormalization factor. 
The black square symbolizes the counterterm proportional to £4. Not all possible diagrams of 
this kind are shown. 
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Figure 6: Next-to-leading order tree diagrams for tt 7 — t- tt tt^tt*^. The black square symbolizes 
the chiral-invariant counterterm proportional to ^1,2,4- 



while diagram (IX) leads to constant amplitudes: 
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(20) 



Fig. 5 shows diagrams with self-energy insertions on external or internal pion-lines. Their 
overall effect is to multiply the tree amplitudes Ai ^2'^ with three times the pion wave function 
renormalization factor Z^ — 1 (see eq.(31) in ref. [23]): 
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Note that we have combined the contribution from the pion-loop proportional to ^ + ln(r?27r/A) 
with that from the counterterm proportional to £4 = £4 + ^/Svr^. 

Fig. 6 shows the remaining tree diagrams for 7r~7 — )■ 7r~7r°7r° generated by the chiral coun- 
terterm Lagrangian C^^l. We use the form of C^^l as written in eq.(20) of ref. [23]. Diagram 
(XI) with a vertex correction at the pion-photon coupling gives rise to the amplitudes: 
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which are obviously proportional to the tree amplitudes. Diagram (XII) involves the four-pion 
vertex from C^^l- For our choice of the interpolating pion field it has no term proportional to 
the low-energy constant £3. With that simplification one finds from diagram (XII): 
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Finally, there is diagram (XIII) involving the (next-to-leading order) 747r contact vertex from 
C[^l- Its polynomial contribution to the amplitudes: 
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is exceptional in the sense that Ai and A2 are not equal. Up to that point the compilation of 
the production amplitudes Ai^2 for the process 7r~7 — )■ 7r~7r°7r° as they arise from chiral loops 
and counterterms is completed. There is still one issue which has to addressed, namely the 
renormalization of the squared pion decay constant fl and squared pion mass ml by chiral 
loops and counterterms (see herefore e.g. eqs. (29,30) in ref.[23]). The chiral Lagrangian £,^77 
operates with their leading-order values whereas the tree amplitudes in eqs. (3,5) have been 
written already in terms their physical values. In order to correct for this difference of order 
(m7r/27r/7r)^ one has to add an extra renormalization contribution of the form: 
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A first crucial check of our calculation is provided by the fact that the ultraviolet divergence ^ 
drops out in the total sums for Ai and A2. Next, one can further simplify the expressions for 
the amplitudes by introducing via the relation: 

^' ^^ +2\n^), 7i = i, 72 = 1, 13 = —, 74 = 2, (27) 
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the "barred" low-energy constants ij which subsume the chiral logarithm ln(m^/A). After 
summing up the expressions in eqs. (18-26) and adding the terms proportional to the chiral 
logarithm ln(m^/A) in the loop amplitudes eqs. (6, 10, 11, 12, 15, 16, 17) one gets the following 
modified and complete counterterm contributions: 
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(477/^)^ 3m^ - s - ti - t2 I 3 3 L 

+ti - 2ss2 + (s + 2si - 2s2 - t2)ti + m^(s - 6si + t2 - 2ti + 6m^ 

-^ml + 2hml{s + 2ml-si-S2)Y (29) 

The remaining finite parts of tlie loop amplitudes are tlien given by eqs. (6, 10, 11, 12, 15, 16, 17) 
witli tlie ^ + ln(m^/A) terms deleted altogether. 

Let us also mention that the other three terms in the chiral Lagrangian C^^l proportional 
to the low-energy constants ^5,6,7 [13] do not contribute to the 7r~7 — ?■ 37r processes considered 
in this work. The £7 term breaks isospin symmetry and is ignored for this reason. The /s 
term requires (for electromagnetic process) at least two external photons. The ^g term gives 
rise to a correction to the pion-photon coupling, which however vanishes for real photons 
ik"^ = = e ■ /c). Let us remind that (i^ is related to the mean square charge radius of the pion 
via: (r^) = (Iq — l)/(47r/7r)^ + Oim^) [24]. The associated 747r vertex vanishes also since the 
commutator in ti{T^[d^U,di,U'^]) terminates at the quadratic order in the pion-field (for our 
choice of interpolating pion-field). 

2.2 Results for 7r~7 — )► tt'tt^tt^ 

In this section we present and discuss the results for the neutral pion-pair production process 
7r~7 — 7> -K^-K^-K^ at next-to-leading order in chiral perturbation theory. We use for the low- 
energy constants Ij the values: li = -0.4 ±0.6, 4 = 4.3 ±0.1, £3 = 2.9 ±2.4, I4 = 4.4 ±0.2, as 
determined (with improved empirical input) in ref. [1]. Applying the flux and symmetry factors 
the total cross section o"tot('S) is obtained by integrating the squared (transversal) T-matrix over 
the three-pion phase space: 

atot(s) = ^^ n,,, ^ duidu2 / dx d(f) kx {qiAi ±^2^2) • (30) 

32rc^f^{s-mi) J J J-i Jo 



Here, uJi^2 are the center-of-mass energies of the out-going neutral pions and qi^2 = y 1^1,2 ^ '"^^• 
In terms of the directional cosines x, y, z the squared cross products in eq.(30) take the form: 

{k X qif = ql{l - x^) , {k X q2f = ql{l - y^) , {k X Qi) ■ {k X ^2) = q^z - xy) , (31) 

together with the relations: 



2 
The Mandelstam variables Si, S2, ti, t2 follow as: 



qiq2Z = U1U2- \/s{uJi + uj2) -\ ^~^ ' y = xz + ^{1 - x'^){l ~ z"^) cos(f) . (32) 



Si = s + ml - 2y/sU2 , S2 = s ± m^ - 2y/suji , 

2 2 

ti=m^H 2^^— (cji -gix) , t2 = m^ H ^^-^ {uj2 - q2y) ■ (33) 



The more definite expression of the tree- level cross section given in eq.(2) is very helpful for 
checking the numerical accuracy of the four- dimensional integration involved in eq.(30). 

Fig. 7 show the total cross section (Ttot(s) for the reaction 7r^7 — ?> tt^tt^tt^ in the low-energy 
region from threshold ^/s = 3171,, up to ^i = 7171,,. The dashed line corresponds to the tree 
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Figure 7: Total cross section for the reaction n^'-f — )■ tt^tt^tt^ as a function of the center-of-mass 
energy y/s. 

approximation and the full line includes in addition the next-to-leading order corrections from 
chiral loops and counterterms. The dotted line follows if the charged and neutral pion mass 
are distinguished (134.977 MeV = -m^o < m.^ = 139.570 MeV) in the tree-level amplitude 
and the three-pion phase space integral (see eq.(2)). As it could be expected such isospin- 
breaking effects are very small. By inspection of Fig. 7 one observes that the total cross section 
for 7r~7 — )■ 7r~7r°7r° gets enhanced sizeably (by a factor of 1.5 — 1.8) after inclusion of the 
next-to-leading order chiral corrections. Although the dynamics of the whole process is much 
richer this feature can be understood (in an approximate way) from the tt+tt^ — ;> tt^tt*^ final 
state interaction. The tt+tt^ — )> tt^tt" interaction strength at threshold is determined by the 
difference of the isospin-zero and isospin-two S-wave vrTr-scattering lengths. When considering 
the corresponding one- loop expression [24]: 



3m7r 
327/J 



1 



m: 



sevrV. 



2 /2 



8 2 8 



(34) 



one sees that the correction to 1 inside the square bracket amounts to about 0.20 (inserting 
the central values of ij). Indeed, the square 1.2^ = 1.44 is close to the enhancement factor 1.5 
of the total cross section at ^i = ArriT,. 

It is also important to investigate the uncertainties which are induced by the present er- 
rorbars 6ij of the low-energy constants ij. Taking the total cross section at ^/s = 6m^ as a 
measure one finds relative uncertainties of: ±5.1% from Sii, ±1.4% from Si2, ±0.2% from Si^, 
and ±0.9% from ^£4. It is comforting that the badly known low-energy constant £3 = 2.9 ±2.4 
has very little infiuence on the observables considered here. The largest uncertainty is actually 
connected with ii and adding them all in quadrature one gets a total relative uncertainty of 
±5.4%. This amounts to a fairly accurate prediction. 

In addition to the total cross section there are more exclusive observables such as the two- 
pion mass spectra. The tt^tc^ invariant mass fi is defined by /i^ = {qi + ^2)'^ = s — Si — S2 + 3m^ 
and it varies over the kinematically allowed range 2171,^ < fi < ^/s — m^^. In tree approximation 
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Figure 8: 7r°7r° mass spectrum for the reaction 7r~7 — )■ vr^Tr'^Tr^ as a function of the tt^tt^ 
invariant mass fi. The curves correspond to center-of-mass energies ^/s = (4, 5, 6, 7)mT^ in 
ascending order. 



the 71 71 mass spectrum can be given in closed analytical form: 



da 
djj, 



aw/i 



Ami 



16vr2/.^( 



m. 



2^3 



{^^' 



m. 



,2\2 



s + m. 



2, , s + m 



/i^ 



W 



2m, 



(35) 



with the abbreviation W = [s — {fi + itLj^)'^] [s — (yU — m,^^)^]. In order to obtain dcr/dfi in general 
one introduces new energy variables w+ = Wi + U2 and a;_ = (t^i — (jJ2)/2 such that /i^ = 
2y^w_i_ + m^ — s. The (ia;+ integration in eq.(30) is omitted and a normalization factor ^ij \fs 



is applied. The lower and upper limit for the duj^ integration are =F\/(Ai^ — 4m^)14^/(4/iy^). 
Fig. 8 shows the calculated 7r°7r° mass spectrum m^ da/djj,. We have multiplied it by the 
constant factor m^ in order to keep the units (/ib) of a cross section. The four pairs of (full 
and dashed) curves correspond to sections at center-of-mass-energies ^/s = (4, 5, 6, 7)m,T^ in 
ascending order. In essence Fig. 8 reproduces the enhancement of the total cross section by 
the next-to-leading order chiral corrections. No further specific dynamical details to which the 
^0^0 jj^g^gg spectrum could be selectively sensitive are visible. 

The alternative combination is to couple one of neutral pions with the out-going tt^. The 
71^71^ invariant mass is ^/si (or equivalently y^). Making use of the relation si = s + 
m^ — 2y^W2 the vr^vr" mass spectrum da /d^fsi is obtained by omitting the duj2 integration 

in eq.(30) and applying an additional factor Jsi/s. Fig. 9 shows the calculated tt^tt" mass 
spectrum m,.,^ da jd^fsl. Again, it is only the enhancement of the total cross section which can 
be inferred from the comparison of the full and dashed curves in Fig. 9. The shape of the vr^vr" 
mass spectrum (i.e. its dependence on ^Js{^ does not distinguish the tree-approximation from 
the full calculation in a noticeable way. 
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Figure 9: tt^tt^ mass spectrum for the reaction n^'-f — )■ 7r^7r°7r° as a function of the tt^tt" 
invariant mass y^. The curves correspond to center-of-mass energies y/s = (4, 5, 6, 7)m^ in 
ascending order. 

3 Charged pion-pair production 

In this section we perform the same calculation and analysis for the charged pion-pair produc- 
tion process: 7r~(pi) -|- 7(A;, e ) — )■ 7r+(p2) + 7r"(gi) + 7r^(g2)- By assigning the four-momentum 
P2 to the out-going positively charged pion 7T^{p2) we can exploit the complete equivalence 
to the 7r~7 — )■ 7i~7i^7i^ reaction concerning its kinematical description. In Coulomb-gauge 
{e ■ pi = e ■ k = 0), eq.(3) constitutes the general form of the T- matrix for 7r~7 — )> tt+tt^tt^ 
and the corresponding Mandelstam variables are defined as in eq.(4). The interchange of the 
two identical vr" in the final state is now described by (si ^ S2, ti -H- 12). The three non- 
vanishing tree diagrams for 7r~7 — )■ ti^ti~ti~ are shown in Fig. 1 and their evaluation leads to 
the following tree amplitudes: 



A\ 



(tree) 



s + m; - Si - S2 s 



si - S2 + t2 



3m?. — s — ti — t2 



U 



mt 



1, 



(36) 



^(tree) 



3m?. 



m^ - Si - S2 s 



Sl - S2 + tl 



U -U 



U 



mt 



(37) 



One can see from the denominators of Ai and A2 ^^° that two diagrams contribute to each 
amplitude. 

3.1 Amplitudes from chiral loops and counterterms 

Beyond leading order the dynamical content of charged pion-pair production 7r~7 — ?> tt+tt^tt^ 
is considerably more extensive than that of neutral pion-pair production 7r~7 — > 7r^7r*^7r° be- 
cause the photon can now couple to all three out-going (charged) pions. Many more diagrams 
with chiral pion-loops and counterterms do contribute. We have evaluated them individually 
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Figure 10: One-pion loop diagrams for vr 7 — ?■ tt+tt tt with three possible couplings of the 
external photon. 



(lib) 






Figure 11: One-pion loop diagrams for tt 7 — )■ tt+tt tt with three possible couplings of the 
external photon. 

and checked the exact cancellation of ultraviolet divergences ^ in the total sums for the ampli- 
tudes Ai and A2. Without loss of information we can restrict the presentation of the analytical 
results to the finite parts of the pion-loop diagrams and to the complete counterterm contribu- 
tion (reexpressed in terms of the low-energy constants ij which subsume the chiral logarithm 
In(m7r/A) generated by the pion-loops). 

The three possible couplings of the external photon for loop diagram (I) (see Fig. 2) are 
distinguished by labelling them (la), (lb), (Ic) and are shown in Fig. 10. Omitting the terms 
proportional to ^ + In(r7i^/A) the corresponding finite parts read: 



A 



(la) 



£^^ 



s + mz 



Sl - S2 



X J{3ml + s - Sl - S2) 



(2s - 2ml -Si-S2 + ti + ts) 

1 
" 2 



(3^ 



A\ 



(lb) 



(4vr/0 



■Sl — S2 ~r 02 , 2 \ 
[S -m^- Sl- S2+ti + t2 

ti -mi 



X 



J{ml + s - Sl - S2 + ti + 12) 



A, 



(Ic) 



Af"^ {Si ^ S2, h ^ t2) 



(39) 
(40) 



Note that we do not list vanishing contributions to Ai or A2 from a diagram under con- 
sideration. Fig. 11 shows the one-pion loop diagrams (Ila), (lib), (lie) obtained from the 
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rescattering diagram (II) (see Fig. 2) by attaching the external photon in the three possible 
ways. One finds the following finite parts: 

,2 



^(Ila) ^ ^(Ila) 



1 r Tfi 

— -^(21s + 7si - 20s2 + 28ti + ta) 

— ti — 12 I O 



3(47r/^)2 3ml- s-h- 

+ 7(s -S2 + ti)(5si + 16s2 - 23s - 18ti - Tts) 
o 



Qlml 
6 



+ 



m^(19m^ - 17s - 2si + 15s2 - 19ti - 2^2) + (s - S2 + ti 



X (7s — si — 5s2 + 6ti 



2t2 



J{2ml - s + S2 - ti] 



(41) 



^1- 



(iib) 



1 1 r TTi 2t71 

TTT^-TT^T ^ ^ + ^(6si + 8S2 - 6s - 3ti - 6t2) 

3{A7rf^yti -mil 2 3 

+ — (5s - 5Si - 16S2 - 2ti + 5t2) + S2(5S2 + Si + ti - s - t2) 

+m2 (4s + m^ - 4si - 9s2 + 2ti + 4t2)l "^(^a) j , 



(42) 



A 



(lie) 



1 



1 



3{4rf^)H2-ml{ 3 



m: 



'lis + 7si - IIS2 + llti - 8t2 



25m,^ 1 

-^ + 77(5 - S2 + ti)(5si + I6S2 + 2t2 - 16s - 16ti) 

o o 

+ U^(9s2 + 4t2 - 9s - 2si - 9ti + 7ml) + (s - S2 + h] 



X (5s - si - 5s2 + 5ii - t2) J{2ml - s + S2 - h) 



(43) 



The additional contributions from the loop diagrams (Ilia), (Illb), (IIIc) with crossed out-going 
7r~ lines (see Fig. 2) follow immediately via the substitution qi ^ q2 as: 



A 



(Ilia) _ A (Ilia) 

— /I2 



A- 



(Ila) 



^(Illb) ^ ^(IIc) 



Si ^ S2, tl ^ t2) , 

AS"''^ (Si ^ S2, tl ^ t2) • 



(44) 



[Si f^ S2, tl ^ t2) , Ar^' = Ar' (Si ^ S2, tl ^ t2) . (45) 

Next, we come to the irreducible one-pion loop diagrams (VI), (V), (VI) with internal photon 
coupling shown in Fig. 3 and interpreted now as diagrams for 7r~7 — ?■ 7r"'"7r~7r~. One finds the 
following finite parts: 



^(iv) _ ^(iv) 



2(s-Si-S2)+ti+t2 (1 



(47rA 



1 



2 2ml - tl - t2 



[s - m^ - Si - S2 + tl + t2) J{ml + S - Si - S2 + tl + t2 



+ (si + S2 - s - ml) J{3ml + s - Si - S2) 



(46) 



A\ 



(V) 



1 



3(47rA) 



13s 
10ml -3s- 3t2 + — r^ + 6m!( 1 



Sl 



G{2ml - s + si - t2) 

3(s-si+t2-2m^)^ 

s + t2 - 2ml 



Gis, 
,2 



6 ' -■■^^V' s + t2-2m2 

(3s - 2si + 3t2 - lOml) J(si) 

J(si) 



J(2m^ - s + si - t2 
15 



(47) 



A 



(V) 



1 



3(47rA)n6 



(19m^ + 13si + 9ti - 2t2 - 2s) + 



3si(2m2-ti-t2) 



+6m: 
+ 



m: 



ti 



+ 



Si(ti 



s + t2 - 2m2 (s + t2 - 2m2) 



si + 2m2 + 



3s?(s-ti 



+ 



3si(ti 



^s + t2 - 2m2)2 s + t2- 2ml 
s — 2si + t2 — 6m^) J(2m^ — s + Si — t 



2(s + t2 - 2ml) 

G{2ml -s + si-t2)- G{sj 

s) 

J{2mi -s + si-t2) 



A\ 



(VI) 



A. 



(V) 



[si ^ S2, h ^ t2 



A 



(VI) 



A\ 



(V) 



>! f^ S2, ^1 ^ ^2 



(4J 



(49) 



where the contributions from diagram (VI) are obtained from those of diagram (V) by applying 
the crossing transformation qi -^ q2- The set of next-to- leading order corrections to 7r~7 — )■ 
TT^TT^TT^ is completed by the total counterterms contribution which reads: 



A\ 



(ct) 



1 



s + 3m_ - 2s2 + ti - ^2 - 



is- Si +t2) 



[s2 — 2m, 



2n2 



(47r/.)2 i 3 
(s - si + t2Y + {S2 + t2 - 3ml)^ 



tl 



mt 



3ml - s -ti-t2 
4s2(s + ml - Si + t2 



3sl 



Ami 



tl -ml 
2{s - si)(2s2 + 2t2 - 5m^) - 3(s -Si + t2- mtV + t 



5s + 5m^ — 4si — 6s2 + 3ti — t2 
+ 



3(s-si + t2)2 1 



+S2(10m2 - 3s2 - 2t2 



+ 



ml 



2\2 
4 



3m2 — s 
-8ml 



tl-t2 



+ 



m. 



+2ml£i 



s + mt 



Si - S2 



tl — m^ ■ 3m2 

S - Si - S2 + t2 



_3mti — s 



^(ct) 



tl — t2 tl 

(ct) 



mt 




Af (Sl ^ S2, h ^ t2) . 



(50) 



(51) 



The relation A2 = ^i|(si ^ S2, ti ^ t2) holds also for the total sum of the loop amplitudes 
and it applies to both reactions 7r~7 — )■ tt+tt^tt^ and 7r^7 — > tt^tt^tt^ in the same way. 

3.2 Results for tt 7 — > tt+tt tt^ 

We are now in the position to present numerical results for the charged pion-pair production 
process 7r~7 — )■ 7r+7r~7r~ at next-to- leading order in chiral perturbation theory. The formula 
for calculating the total cross section o"tot(s) is given in unchanged form by eq.(30). We use 
consistently the same values: li = -0.4 ± 0.6, I2 = 4.3 ± 0.1, I3 = 2.9 ± 2.4, £4 = 4.4 ± 0.2, for 
the low-energy constants ij as in subsection 2.2. 

Fig. 12 shows the total cross section for o'tot(s) for the reaction 7r~7 — > tt+tt^tt^ in the low- 
energy region from threshold ^/s = 3m^ up to y/s = Tm,,^. The dashed line corresponds to the 
tree approximation and the full line includes in addition the next-to-leading order corrections 
from chiral loops and counterterms. By inspection of Fig. 12 one observes that the total cross 
section for 7r^7 — )> tt+tt^tt^ remains almost unchanged in the region ^/s < Qm.,^ after inclusion 
of the next-to-leading order chiral corrections. This striking result is in marked contrast to the 
behavior of the total cross section for neutral pion-pair production 7r~7 — ?> -k^-k^-k^ (see Fig. 7) 
Although the dynamics of the whole process is much richer this feature can be understood (in 
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Figure 12: Total cross section for the reaction 7r^7 — > tt^tt^tt^ as a function of the center-of- 
mass energy ^/s. 

a suggestive way) from the tt^tt^ — )■ tt^tt^ final state interaction. By considering the one-loop 
expression for the isospin-two S-wave tttt scattering length [24]: 



m^ 



167r/2 



m: 



12vr2/2 



'i + 24 - 



3£, 



3^4 3 
^+8 



(52) 



one sees that the correction to 1 inside the square bracket amounts to the very small number 
—0.017 (inserting the central values of £j). Chiral corrections (even at two-loop order [1]) 
affect the isospin-two vrvr-interaction only very weakly and this feature seems to be reflected 
by (Ttot(s) in Fig. 12. Note however, that the argument made here is only suggestive and not 
rigorous, because the on-shell tt^tt^ — )> tt^tt^ final state interaction does not factor out of 
the production amplitudes Ai and A2 in an obvious way. The same caveat applies to the 
7r"'"7r~ — )■ Tr'^Tr'^ final state interaction which has been used as an argument for the observed 
enhancement in subsection 2.2. 

Let us also comment on the theoretical uncertainties which are induced by the present 
errorbars 6ij of the low-energy constants ij. Taking again the total cross section at ^/s = 6112^^ 
as a measure one finds relative uncertainties of: ±4.8% from 6ii, ±1.6% from S£2, ±0.3% from 
(5^3, and ±1.0% from (5^4. As for the reaction 7r~7 — > 7r~7r°7r° the largest uncertainty goes along 
with ii and adding them in quadrature one gets a total relative uncertainty of ±5.2%. This 
amounts again to a fairly accurate prediction. 

The more exclusive observables than the total cross section are the two-pion mass spectra. 
Fig. 13 shows the calculated tt^tt" mass spectrum 171^^ da/dfi as a function of the vr'vr" invariant 
mass fi. The four pairs of (full and dashed) curves correspond to sections at center-of-mass- 
energies y/s = (4,5,6,7)m7r in ascending order. In essence Fig. 13 reproduces the features of 
the total cross section, namely a slight enhancement above ^/s = Gm^r by the next-to-leading 
order chiral corrections. The tt+tt^ mass spectrum m^ da/d^/sl shown in Fig. 14 indicates 
some more interesting structures. The dip of the tt+tt^ mass spectrum at intermediate tt+tt" 
invariant masses ^/Sl becomes much more pronounced when including the next-to-leading order 
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invariant mass /x. The curves correspond to center-of-mass energies ^Js = (4, 5, 6, 7)111,, in 
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chiral corrections. This distinctive feature holds e.g. at the center-of-mass energy ^/s = 5m„ 
where the total cross sections in tree and one-loop approximation are equal. The tt^tt^ mass 
spectrum of the reaction 7r~7 — )■ tt+tt^tt^ therefore seems to provide an interesting indicator for 
the role of chiral (pion-loop) dynamics beyond leading order. It is expected that the upcoming 
high-statistics data of the COMPASS experiment at CERN can reveal such dynamical details. 
Of course, the squared (transversal) T-matrix \k x {qiAi + 5*2^2) P with its full dependence on 
pion energies and angles incorporates still much more dynamical information. 

In passing we note that ^/s = (6 — 7)m„ is presumably the maximal center-of-mass energy 
up to which a one-loop calculation of the processes 7r^7 — )> Svr in chiral perturbation theory can 
be trusted. At still higher energies the contributions from meson resonances (such as ai(1260), 
02(1320) etc.) will start to play a prominent role. In the context of such considerations it 
should be kept in mind that the effects of the resonance tails at low-energies are encoded in 
the empirical values of low-energy constants ij. The role of the low- lying p(770) resonance 
occurring in the isospin-one 27r-subsystem needs to be investigated by studying an appropriate 
resonance model for 7r^7 — )■ Svr. Respecting fully gauge-invariance in the construction of such 
a resonance model (with inclusion of finite resonance widths) represents some challenge. 
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